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Abstract 

We present a study of the magnetic order and the structural stabihty of two-dimensional quantum 
spin systems in the presence of spin-lattice coupling. For a square lattice it is shown that the 
plaquette formation is the most favourable form of static two-dimensional dimerization. We also 
demonstrate that such distortions may coexist with long range magnetic order, in contrast to the 
one-dimensional case. Similarly, the coupling to Einstein phonons is found to reduce, but not to 
eliminate the staggered magnetic moment. 

In addition, we consider the renormalization of the square lattice phonon spectrum due to spin- 
phonon coupling in the adiabatic approximation. Towards low temperatures significant softening 
mainly of zone boundary phonons is found, especially around the (vr, 0) point of the Brillouin zone. 
This result is compatible with the tendency to plaquette formation in the static limit. We also 
point out the importance of a "magnetic pressure" on the lattice due to spin-phonon coupling. At 
low temperatures, this results in a tendency towards shear instabilities of the lattice. 

PACS numbers: 03.65.-w, 73.43.Nq, 03.75.Lm, 32.80.Bx, 05.70.Fh 
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I. INTRODUCTION 



Considering the immense wealth of materials newly synthesized or found in nature with 
low dimensional magnetic structure, one finds that the simple Heisenberg Hamiltonian, 
widely accepted as the paradigm of quantum magnets , is often complicated by further 
interactions. Many compounds thus show ground states and phases, different from the 
generic long range magnetic order of Heisenberg systems. 

In this paper we study the coupling of the magnetic degrees of freedom to lattice vibra- 
tions as a prominent example of an interaction "beyond the Heisenberg-Hamiltonian" . This 
coupling, often referred to as " spin-phonon" interaction is generic to all magnetic materials 
and thus the problem to what extent the magnetic structure is changed by its presence 
concerns all low dimensional magnetic compounds. 

Though spin phonon couplings are ubiquitous in nature, the possibilities of treating them 
with some theoretical rigor are limited. The main reason is that magnetic and phononic 
excitation energies are not well separated. In contrast, for electron phonon coupling we 
can dwell on the fact that the bulk of the electronic excitations lies much higher in energy 
than the phonon modes. This allows, in general, to use the adiabatic approximation when 
calculating the renormalization of phonon frequencies due to electron-phonon coupling. For 
spin-phonon coupling, however, this concept cannot be applied, in general. 

One approach quite often used is to include local, Einstein-type phonons as quantum 
mechanical objects. This model may be applied for certain vibrations of ligand atoms around 
the magnetic ions, but it is not very realistic as it violates the infinitesimal translational 
invariance of the lattice. 

For relatively low-lying acoustic phonons, which are dominated by the heavy magnetic 
ions, the response of the spin system may be treated in adiabatic approximation. In one 
dimension the paper of Cross and Fisher 3| has treated the spin-phonon coupling along such 
lines. In higher dimensions no study of spin-phonon coupling in the adiabatic limit has so 
far been presented. 

In addition any effect of the magnetic pressure on the phonon vibrations has been ignored 
so far. The magnetic energy of a square lattice of spins with INN antiferromagnetic coupling 
J {a) is given as = U{T)J{a) which is typically of the order —0.1 eV per atom. We may 
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compare this magnetic energy with the energy EMad = —C(m^ of a cubic lattice of positive 
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and negative elementary charges where the Madelung constant aM is of order of unity. For 
rock salt we find e.g. Ejuad ~ —8eV. However, the pressure |^ for the electrostatic system 
is ~ 0.1e\^/A^, while ~ 0.015e\^/A^ for the magnetic system. Although there is a factor of 
~ 100 difference in the binding energies, the magnetic pressure is only a factor 10 smaller 
compared to the pressure of the electrostatic system. The main reason for the relatively 
large magnetic pressure is the large spin phonon coupling caused by the strong distance 
dependence of the super exchange J{a), which cannot be neglected in describing realistic 
quantum spin systems. 

As an important special case we also discuss the limit of very strong spin phonon coupling, 
which leads to a dimerization of the lattice. Typically, the response of the ground state 
energy of a two-dimensional magnetic system to dimerization is of second order and thus 
much weaker than in the one-dimensional case. Two dimensional models with dimerization 
are thus essentially different compared to their one-dimensional analogues, in particular in 
two dimensions dimerization does not necessarily lead to a breakdown of magnetic long 
range order. Also the dimerization pattern, which leads to the lowest ground state is not 
clear in two dimensions. 

This paper is organized as follows: In Sect IIA we study statically dimerized models, 
in particular we determine the optimal dimerization pattern in two dimensions. Sect IIB 
is devoted to Einstein phonons coupled to two-dimensional spin systems. In Sect III we 
calculate explicitly the phonon spectrum and study structural instabilities driven by the 
spin-phonon coupling including the effect of the " magnetic pressure" . 



II. STATICALLY DIMERIZED MODELS AND EINSTEIN PHONONS 
A. The dimerized model with optimal deformation patterns 

To approach the problem we assume in the first place that in two dimensions similar 
to one dimensional systems the spin lattice interactions lead to a transition to a dimerized 
state i.e. to a pattern of strong and weak bonds and we address the problem of the two 
dimensional correspondence of a dimerized chain. 



Stair, plaquette and meander configurations (see FigH]) are o 
els, which have been discussed controversially in the literature 



)vious choices for such mod- 
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FIG. 1: Plaquette, stair, and meander configurations as possible dimerization patterns with mini- 
mal unit cell in two dimensions. 

Here we show by a straightforward Monte-Carlo calculation [3] that the plaquette-models 
have lowest energy. This is convincingly demonstrated in Fig|2l where the extrapolated 
energies of the three configurations are shown. 

It is not clear however, whether configurations with larger unit cells or even disordered 
systems can still have lower energies. To clarify this point, we analyze Hamiltonians of the 
type 



Bij) = 0. Here Aij and Bij specify the pattern and 6 gives the strength of the deformation. 
This Hamiltonian of course includes the configurations of Fig.l as special cases. 

This more general study stems from an expansion of the free energy of the system Eq. ^ 
in 6 



Writing down a(T) explicitly we observe that it can be viewed as Hamiltonian of a two layer 
Ising model with Ising spins Ay and B^j 




for all possible Aij, Bij = ±1, with equal number of strong and weak bonds M = Xli + 



F{6, T) = F(0, T) + \a{T)5^ + y^KT)5' + 0{5'). 



(2) 




(3) 



+Kyy{di,d2) B,,B,+d„j+d, 

i+di,j+d2 

+Ky''{di, d2)BijAi+di,j+d2} 



where the couplings of the Ising model are given by dynamic dimer-correlations 




(4) 
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FIG. 2: Extrapolated ground state energies for the three dimerization patterns shown in Fig^ 



evaluated in the two-dimensional Heisenberg model, i.e. for 5 = 0. Here q and r are either 
X or y corresponding to dimer operators D^j = SijSi+ij or Dfj = SijSij+i. This means 
the quantum nature of the model is incorporated in the long ranged Ising couplings which 
depend on the Euclidean dynamical dimer correlation functions of the isotropic Heisenberg 
model. This "perturbative approach" is somewhat restrictive but we are interested in the 
phenomenologically relevant small dimerizations. Also, with some more numerical effort 
quadruple and higher correlations could be studied and for large dimerizations the prob- 
lem becomes rather trivial, since it is reduced to one dimensional Heisenberg chains whose 
properties are well known. 

The dimer correlations Eq. 0] of the Heisenberg model (as well as the data shown in 
FigEI) were evaluated using a quantum Monte Carlo fQMC) loop algorithm 0,0], based on 
a path integral representation of the partition function [l(j|. Applying the loop algorithm, 
autocorrelation effects play a minor role due to global spin updates, and expectation values 
of diagonal and off- diagona l operators are calculated efficiently within the framework of 
i_d e*n,ato.s QQ. Furthe, by t^.„g the continuous tin,e version of the ^gortthur 
|12| . finite size effects in Trotter direction are avoided. 

With the dimer correlations Eq|3] as coupling constants of the Ising model we performed 
a standard classical Monte-Carlo simulation together with some cooling procedure QjQ- 
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This analysis gives clear evidence for crossing stripe patterns in the plane of the A and B 
spins as ground state of the Ising model This amounts to a plaquette structure as 

lowest pattern of the model Eq^J 

In contrast to the one-dimensional case where a{T) is divergent when T approaches 
zero, a(T) stays finite in two dimensions, which a posteriori justifies our ansatz. 

A crucial point which further distinguishes the one dimensional case from the two- 
dimensional is the presence of both long range magnetic order and a finite dimerization 
in the two-dimensional plaquette-model, whereas in one dimension any finite 6 leads to a 
spin liquid ground state with a gapped excitation spectrum. 

For static phonons we definitely observe a breakdown of long range magnetic order for 
sufficiently strong deformation 6. The critical value 6c takes rather large and in general 
different values for different patterns. For the stair configuration 6c takes its maximal value 
1, at which the system decouples into one-dimensional subsystems. For the plaquette config- 
uration, 6c is considerably smaller than 1 (in the cluster series expansion approach 6c ~ 0.3 
jlS^). while we determined from our data 6c = 0.291(3). 

B. The effect of Einstein phonons in two dimensions 

To study the impact of a spin-phonon coupling on the magnetic long range order we 
next consider a two dimensional model coupled to Einstein phonons. Although the model 
is not realistic in many respects, it has still attracted a lot of attention and for the physics 
of its one dimensional form a fairly clear picture has emerged by now. In particular it 
is well known, that the quasi long-range order, leading to a logarithmic divergence in the 
structure factor at g = vr is destroyed by a relatively small spin-phonon coupling. This rises 
the problem to what extent the strong long range order of 2d systems is influenced in the 
presence of Einstein phonons. To answer this, we calculate the expectation value of the 
staggered magnetization operator, 

x,y 

which is a measure of Neel-order in the ground state. We adopt the common procedure (see 
{3,13) compute the expectation value of Mfj 

:= ^(OIM^IO), (6) 
6 



0.16 



0.12 




0.08 



0.04 



A g=0. 1 J, 0}=8J: M=0.298(5) 
□ g=(a=0.67J: M=0.275(5) 
• Sandvik: M=0.3070(3) 



0.08 
1/N 



0.16 



FIG. 3: Staggered magnetization for the model Eq[7| as function of linear lattice size. The 
black circle is the value for the Heisenberg-model from Ref. 2ol |. 



for the full quantum-Hamiltonian of spins Sij at sites ij coupled to Einstein phonons b, 
and b^j by a spin phonon coupling g 



(7) 
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221, we find, as 



Employing the Monte-Carlo method, which was developed in Ref. 
expected, Heisenberg like behaviour in the case of small spin-phonon couplings g and large 
values of the phonon frequency u. It is remarkable however that for parameters for which 
in the one-dimensional model one finds clear evidence for a finite correlation length at T=0, 
the staggered magnetization is again nonzero, and its extrapolated value is only about 10% 
reduced compared to the two-dimensional Heisenberg model. Also, displays a dominant 
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24| derived from the nonlinear sigma model description 



1/A^- finite size behaviour 
of the Heisenberg model. 

Though we cannot exclude a breakdown of long range order by studying even larger 
spin-phonon couplings or different, maybe more elaborate coupling mechanisms, we find, 
that a realistic coupling strength comparable to the ones found in one dimensional magnetic 
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materials will not lead to a breakdown of long range order. So, the new feature we expect 
in two dimensions is the possibility of a structural phase transition driven by spin-phonon 
interactions without breakdown of long range order. 

III. STRUCTURAL PHASE TRANSITIONS 

To investigate the structural phase transition also from a phenomenological point of 
view, we consider the following two-dimensional Heisenberg-Hamiltonian coupled to classical 
phonons 

H = '^JY1 l[^ + Kui-Um)R^SlSm, (8) 
l,m 

with Si a spin 1/2 operator at position ri, Rim = Ri — Rm the distance vector between 
sites / and m at equilibrium. Rim = Rim/\Ri — Rm\ and ui = fi + Ri the displacement vector 
of the (l)th spin from equilibrium. The summation in Eq. |H1 runs over nearest neighbors 
1 and m. (Note that to simplify the notation we abbreviate from now on the two spatial 
indices used e.g. in Eq. ^by only one index and continue to use the coupling strength 
2 J per bond.) 

The total Hamiltonian consists of a sum of the Hamiltonian H for all layers plus the 
phonon contributions $, which we assume to be derivable from a simple model with central 
force potentials $i(|r« — ^ml) and $2(|n — "^ml), which depend on the nearest neighbor and 
the next nearest neighbor (diagonal) distances between the ions 

^ = ^ E ("^id^"' - ^"-D + "^^dn - rm\)) . (9) 
Thus in equilibrium the total energy E{a) per site is given by 

E{a) = 2$i(a) + 2$2(v^a) + 2U{T)J{a) (10) 

where U(T) < is the internal energy per site of the Heisenberg model with J=l and a 
is the lattice constant. The magnetic energy U{T)J{a) has the tendency to compress the 
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lattice as the super exchange couphng J (a), which may be derived from a Hubbard- type 
model, increases with decreasing lattice constant a. When we treat this magnetic pressure 
as an external pressure acting on the lattice, the equilibrium condition ^ = yields 



2B + 4B' + -JXU{T) =0. (11) 



Here A = '^^^"^ is again the spin phonon coupling, and following [25[ the force constants 
B and B' are defined as 



1 

:= (12) 

Rim dRlm 



with Rim = \Ri — Rm\ the equilibrium separation for either first or second neighbor pairs. 

(i) 



The full force constant matrix $1^*1 is given by 



$S = ^^"^RtmRL + B'^\6'^^ - RtmRL) (13) 

with 



^^^^ = ;7^*^^n^^n.)- (14) 

dRlm 

In any reasonable model of INN and 2NN force constants A, B and A\ B' 

\B\ < A and \B'\ < A' (15) 

can be expected from the requirement that longitudinal phonon frequencies are in general 
considerably larger than transverse ones. One should keep in mind that our model for a 
square lattice of magnetic ions ignores additional forces acting via the closed shell ligands, 
which may also contribute to the equilibrium condition EqlTT| so that without the magnetic 
pressure, B + 2B' could be chosen positive. When, however, the magnetic pressure is turned 
on, by lowering T, B + 2B' will become negative (or strongly reduced) which will affect, in 
particular, the stability of the lattice against shear. 

The elements of the dynamical matrix GQ^{q) of the "bare" phonons are given by 
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FIG. 4: Eigenspectra Lof ^ of the matrix g for T/J = 0.1, T/J = 0.9 and T/J = 2.0 along a 
triangular path in g-space. 

G^''{q) = 2A{1 - cos x) + 2B{l- cosy) + 2{A' + B'){1- cos X cosy) (16) 
Gl^{q) = 2A{1- cosy) + 2B{1- cos x) + 2{A' + B'){1- cos X cosy) 
Giy{q) = 2{A' - B')smxsmy. 

with X = Qxtt and y = qya. 

We find the following phonon dispersion relations for longitudinal and transverse modes. 
For q = (x, 0) 

m toliq) = 2{A + A' + B'){l-cosx) (17) 
m iuliq) = 2{B + A' + B'){l-cosx). (18) 

for q = (x, x) 

m iuKq) =2{A + B){l-cosx)+4:A'sm'^x (19) 
m uj^q) = 2{A + B){1 - cosx) + AB' sin^ x (20) 

At long wavelength, the transverse frequencies are given as 

mujl{x,Q) = {{B + 2B') + A'-B')]x'^ (21) 
m ujI{x, x) = [2{B + 2B') + A- B)]x^. (22) 
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In addition to the magnetic pressure, which is first order in A, the spin phonon couphng 
renormahzes the "bare" phonon frequencies similar to the renormalization due to electron- 
phonon couphng. 

Similar to the analysis described in the first part of the paper, we again expand the free 
energy of the system EqElin terms of the lattice displacement and find an expression, which 
is similar to EqlHl However here the Aij of EqOl are replaced by mJ+ij — wfj and the Bij 
by u^j+i — u^j- That is compared to Eq. El the classical displacements are no longer Ising 
variables but take continuous values. Thus to second order in the spin-phonon coupling we 
find for the dynamical matrix 

G{q}=Go{q)-X'9{q), (23) 

where the entries of the spin-phonon contribution 

c/^^ = 2(1 - cosx)i^^^(a;,y), (24) 

g^y =2{l-cosy)Kyy{x,y), 

gxy ^ gyx ^ _ ^ix^i^i _ e-iy)K''y{x, y) + c.c. 

involve the Fourier transforms K'^^{x,y) of the dimer-dimer correlations EqEJ Note that 
the matrix g is real and symmetric, and from the symmetries of the dimer correlations 
follows K''y{7i,TT) = 0. 

Typical eigenvalues for the spin-phonon coupling matrix g{x, y) along the symmetry lines 
of the square lattice are shown in FigEl Several features are noteworthy: 

(i) there is a considerable T dependence not only concerning the magnitude of the renor- 
malization, but also the q-dependence. 

(ii) there is a maximum of the renormalization around T/J=l. 

(iii) along (x, 0) there is no renormalization of the transverse branch, similar to the absence 
of the electron-phonon coupling for transverse modes in a nearly free electron model. 

(iv) at small T, a significant splitting of the eigenvalues along (x, x) evolves, indicative of 
an increase of K^y at low T. 

(v) also, the renormalization of the longitudinal branch along (x, 0) shows an increase of the 
effective 3A^A^ force constant towards low T. 
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FIG. 5: Eigenspectra uJi^2 in zeroth (dashed line), first (solid line) and second order (circles) along 
a triangular path in g-space for spin-phonon coupling A = —0.6, A' /A = 0.8, B/J = —0.075, 
B' / J = —0.163 and T/J = 0.1. Note, that in the region of instability with 2 < '^^ plot 

The points (iv) and (v) indicate that the effective forces between the atoms transmitted 
via the spin system become longer ranged at low T. 

The phonon renormalization due to spin-phonon coupling may lead to significant softening 
and even to lattice instabilities. 

The magnetic pressure effect lowers specific transverse phonon modes at long wavelengths, 
i.e. it destabilizes the lattice against certain shear deformations. This can be seen in FigJSl 
where the undisturbed spectrum of a square lattice (zeroth order in A) which follows from 
EqfTBl by setting i? = i?' = is shown together with the spectra including the external 
pressure and the second order contribution. Here the equilibrium condition Eq^^is a strong 
constraint on the system and the shear instability dominates the phonon renormalization. 

We should keep in mind however that realistic planar quantum antiferromagnets usually 
are ternary or even quaternary oxides and that the equilibrium condition of the lattice may 
be dominated by ligand contributions, not included in Eq^2 As a consequence, we may 
ignore the magnetic pressure effect and consider solely the phonon renormalization due to 
the second order effect. Here primarily zone boundary phonons are lowered, as may be 
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FIG. 6: Energies a;i^2 of the matrix G for B = B' = 0. Again for Luf2<0we plot —y/—uJi^2- 

seen from FigEl For the calculation we have have chosen typical values for the model force 
constants, which approximately reproduce the acoustic phonon modes of typical cubic oxides 
with rock salt structure 26| such as MgO or BaO. 

For a wide variety of A' /A ratios we find that the softening is strongest near (vr, 0); only for 
somewhat pathological choices of A, A', B, B' such as A' — B' 0, A' + B' ^ A,we find that 
the (tt, vr) frequency goes unstable first. The (vr, 0) instability would lead to plaquette-like 
distortions, while the (vr, tt) instability would yield stair-case distortions. 

This means that our results from spin-phonon coupling agree with the observation dis- 
cussed above, that the plaquette is the energetically most favorable dimerization pattern. 



IV. SUMMARY 



In this paper we have studied magnetic and structural instabilities of two-dimensional 
quantum antiferromagnets. Of the various dimerization patterns the energetically most 
favourable one at T = is the plaquette order. It was also shown that the various dimeri- 
sation patterns do not necessarily lead to a breakdown of long range magnetic order. In 
contrast to one dimension, the coupling to Einstein phonons is found to reduce, but not to 
destroy the staggered magnetic moment. 
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The tendency to form dimerization patterns at sufficiently large spin phonon coupling 
is also reflected in the large renormalization of speciflc phonon frequencies, in particular 
near (vr, 0) and (tt, tt). Our results indicate that again, plaquette formations are the most 
favourable structural distortions. 

In this paper we also point out to the importance of a "magnetic pressure" on the lattice 
caused by the relatively large spin phonon coupling. This pressure effect builds up with 
decreasing T and may lead to a significant decrease, if not instability of specific shear elastic 
constants. Presently we have treated the magnetic pressure as an external one. In principle 
it can also be included into the calculation of phonon frequency renormalization, e.g. by 
carrying out so called frozen phonon calculations, which is the ultimate form of the adiabatic 
approximation. 

We finally remark that, wherever in the Brillouin zone there occur renormalization effects 
due to spin phonon coupling, also a significant increase of the phonon line width is expected. 
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